Parametric instability induced scalar gravitational waves from a model pulsating 

neutron star 



Charles H.-T. Wang/'^.Q Paolo M. Bonifacio/ Robert Bingham,^' 3.| and J. Tito Mendonga''' =^ 'I 

'SUPA Department of Physics, University of Aberdeen, King's College, Aberdeen AB24 3UE, UK 
^Rutherford Appleton Laboratory, STFC, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
^SUPA Department of Physics, University of Strathclyde, Glasgow G4 ONG, UK 
^ GOLP/Centro de Fisica de Plasmas, Instituto Superior Tecnico, 1049-001 Lisboa, Portugal 

We identify a new dynamical mechanism for a strong scalar gravitational field effect. To illustrate 
this mechanism, we analytically and numerically investigate the parametric excitation and emission 
of scalar gravitational waves by a radially pulsating model neutron star. 
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Introduction — Within the metric description of grav- 
itation, scalar-tensor (ST) theory of gravity provides a 
class of natural extensions to general relativity (GR), by 
including a scalar field as part of gravity. Though a mul- 
tiplet of scalar fields can be considered, here we shall 
focus only on a single real-valued scalar field and denote 
it by (p. It modifies gravity in the following way: In the 
absence of this field, the motion of matter fields is deter- 
mined using the metric g*i^, referred to as the "Einstein 
metric" as per GR. In the presence of the scalar field <j>, 
the motion of matter fields are affected by the "physical 
metric" gab obtained from the Einstein metric through 
a conformal mapping gab = ^{4')'^ Qabx-^^^^^ ^{4') is a 
theory-dependent "coupling function" Since ef- 
fectively re-scales the mass of particles, it is also known 
as the "mass function" [3, [3| ■ The coupling strength be- 
tween matter and scale fields 4> is proportional to the 
field derivative a{(j>) := da{(j))/d4> of the natural loga- 
rithm a{(f>) := In A{(l)) of the coupling function. (See the 
field equations below.) Thepioneering scalar-tensor the- 
ory by Brans and Dicke 0, is equivalent to the simple 
choice of a linear function a{(f>) = ao(j) where ao is a con- 
stant and is related to the Brans-Dicke parameter wbd 
by ckg = (2ci;BD + 3)~^. Thus the larger ujbd, the weaker 
scalar field coupling. The most stringent test to date of 
the Brans-Dicke theory is provided by measuring the fre- 
quency shift of the radio signals to and from the Cassini- 
Huygens spacecraft, suggesting that uj^b > 40000 _^]. 
Further tests could be provided by the measurements of 
(tensor) gravitational wave forms from neutron stars and 
black holes using advanced laser interferometers Jj, ^] . 



lol 11 1 . The implications of this strong- field effect have 
attracted considerable recent attention e.g. [l2, 13, 14 1. 

The purpose of this paper is to point out a new strong 
field effect for /3 > 0. The scenario involves dynamical 
parametric instability for the scalar field inside a rapidly 
pulsating strong source of gravity, such as a young neu- 
tron star. In order to gain some understanding of this 
new effect, an idealized neutron star model is analyzed. 
We show that for some sufficiently large (3 together with 
certain physically reasonable conditions on the star, the 
pulsating model neutron star behaves like an optical cav- 
ity in which resonant scalar waves are parametrically 
amplified. The surface of the star then acts like a par- 
tial (anti-phase) refiector that releases travelling scalar 
waves analogous to how optical laser is generated. This 
work therefore provides an initial estimate of the effect 
that can be extended for further investigation with real- 
istic stars, including the possible energy transfer from a 
collapsed star core to stalled shock waves in supernova 
formations and other astrophysical problems 15 1. 

The scalar gravitational field — The general action for 
scalar-tensor gravity including matter can be expressed 
as dH: 
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where a, 6, • • • = 0, 1, 2, 3, d*x — dx^d^x and x° — ct, and 
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The weakness of the Brans-Dicke type scalar coupling 
can be explained in terms of the relaxation of a more gen- 
eral nonlinear function a{4>) to its local minimum during 
the cosmological evolution The effect amounts to a 
local attractor (j)Q around which a{(j>) « /3((/) — (/)o)^/2 
with some constant (3 > 0. On the other hand, if 0o 
corresponds to a local maximum, then /3 < and "spon- 
taneous scalarization" could occur 0, Indeed, for 
(3 < —4 large deviation of the scalar field from its cos- 
mological value may be developed inside a neutron star 



is the action for the scalar field (j) in terms of a potential 
function V{(j)) and Sm — Sm[fp, A'^{(p)glb\ is the action 
for matter fields in terms of a coupling function A{(j3). 
The Minkowski metric is given by rjab = diag(— 1, 1,1,1). 
The effective stress tensor for 6 follows from (121) as 
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The field equation for (j) follows from varying the total 
action S* in ([T]) as 



dV{<i>) 47rG. 



(4) 



where is the Laplace-Beltrami operator and T, is the 
contracted stress tensor of the matter field with respect 

to 9*ab- 

We shall be concerned with the values of (j) near a lo- 
cal minimum of a{(j)). Thus up to an additive constant, 
equivalent to a re-scaling constant for the metric gati we 
have approximately 



a(0) = i/3(/.2 



(5) 



for some constant /3 > 0. For simplicity we shall consider 
here the quadratic potential 



(6) 



that gives rise to an effective mass mo = fJ-ofi/c of the 
scalar field i/i in vacuum. Then ^ becomes 



— /io0 = U4> 



where 
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The model neutron star — It is interesting to note that 
([7]) is a homogenous wave equation for (p. As such, para- 
metric excitation of (j) through a time varying function 
U could occur. However, generation of large amplitude 
scalar wave does not necessarily follow unless parametric 
instability takes place. Such a condition would require 
rapid variations of a high density gravitational source 
with a sufficiently large value for p. In addition, damp- 
ing of the parametric excitation through radiating scalar 
waves must also be sufhciently small, otherwise the para- 
metric excitation would be stabilized. In this paper, we 
investigate a particular scenario by means of a model 
pulsating neutron star to illustrate how these conditions 
can be satisfied through the parametric instability of the 
quasi- normal modes of the scalar field inside the star. 
We envisage that similar mechanism could be identified 
in other violent events in astrophysics and the early uni- 
verse. 

The radiation of scalar gravitational waves by a ra- 
dially pulsating star using the Brans-Dicke theory was 
studied in [l^. (See also The emission mechanism 
in that case is essentially the same as that for the conven- 
tional tensor gravitational waves but monopole radiation 
is possible. In our present case of quadratic scalar cou- 
pling using ([5]), the dynamical structure is completely 
different, as we will see below. 



Following [l6|, we adopt 



(9) 



as the contracted stress tensor for matter inside a spheri- 
cally symmetric star model with radius R, density p and 
pressure p. Outside the star we have simply p = = p 
and hence = 0. For this model, we shall further as- 
sume c^p ^ p so that ([5]) can be approximated by 



(10) 



For simplicity, we shall also ignore the curvature of g*^ by 
approximating it with the Minkowski metric rjat- There- 
fore the scalar field equation ([7]) reduces: 



(11) 



where A is the 3-dimensional Laplace operator. 

The star is modelled as a solid sphere with a constant 
equilibrium density po- Its density fluctuations are de- 
scribed by the wave equation 



d'oP~{v/cfAp^O 



(12) 



where v is the speed of the density /pressure wave. For 
simplicity we shall consider a single mode radial oscilla- 
tion of the density subject to zero boundary condition at 
the surface r = R so that 

p = po[l- eXm(r) cos(rj„it)] (13) 

for some positive integer m as a mode index, where 

= S- (14) 
K 

is the oscillation frequency, e is a dimensionless amplitude 
parameter, and 



R 

Xn{r) := — sin(K„r) 
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(15) 



with the wave number 



rnr 



for n — 1,2,... The following orthogonality relation 
holds: 
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for any m,n = 1,2,... Substituting (fT5)) into ^TU\i we 
have 



where 



U = Uo[l ~ exm{r) cos(n„t)] 



Uo := 2 — Po - 0- 
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Then pA|) becomes 

dlcj) - A0 + - eUoXrnir) cos(17„,O0 = (19) 
where 

/i2 := ^2 + Uo. (20) 

The energy density and flux of can then be evaluated 
from ^ using g*^ = r]ab to be 



[</',O0,O+??"^'/',a<^,/5 + 2V((^)](21) 

(22) 
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respectively, where a,/? = 1,2,3. The potential takes 
the forms ^(0) = ^^i^tf)^ and ^(0) = |/io'/'^ inside and 
outside the model star respectively. 

Approximate normal modes of the scalar field inside 
the model star — In the absence of the density oscillation, 
i.e. e = 0, the model star has a constant density poi 
inside which cj) satisfies the Klein-Gordon equation with 
mass parameter /i > /iq given in (|20p . For /i//io ^ 1 the 
star surface at r = i? behaves like a perfect (anti-phase) 
reflector for outgoing 0. As such the scalar field inside 
the star can be approximated by a standing wave subject 
to the boundary conditions (t){R^ t) — as follows: 



(23) 



where n — 1,2.... Each 0„ denotes a normal mode with 
^„(i) = 5ft(^„oe-*'^"* (24) 
where ipno is a modal amplitude constant and 
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To evaluate the energy associated with these normal 
modes, we first obtain the energy density of (f> inside the 
model star 



.[(0,o)' + (0,,)2 + AiV^ 



SttG, 



(26) 



by using (PT|) in spherical coordinates. Using 

and (Hil), this yields the following energy for 
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Quasi-normal modes and damping due to transmitted 
scalar waves — For finite /i//io > 1, the star surface does 
allow some scalar wave to propagate across it. 

By taking into account the resulting loss of energy, we 
can refine 0„ in (|23p to be quasi- normal modes. To this 



end, we assume that f24)) is valid over a few circles of 
oscillation at angular frequency uj„. Exterior to the star, 
this yields the scalar field 



hrt r 



(28) 



exterior to the model star {r > R), where 9n is a constant 
phase and 



p -'^^ ,,2 



Using (HOI), (E3 we have 
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(29) 



(30) 



The power carried by outgoing waves then follows from 
as (at r :s> R). 



(31) 



P„ = 47rr2|/'3| = ±-^i?V2^c.„fc„ 
This yields the damping factor dn'- 

, Rn IcP" 

dn = — = 



En RuJnkn 



(32) 



for ifn, which now describes a quasi- normal mode satis- 
fying the damped oscillator equation: 



dVn , , dipn 2 n 



(33) 



Parametric excitation of the normal modes — We now 
estimate the parametric excitation of the normal modes 
of the scalar field in the presence of the density oscillation 
described by (fT3| with e ^ 0. We shall first neglect the 
efi'ect of damping as just discussed and then take this ef- 
fect into account later on. In order to gain a simple rough 
estimate of the excitation effect we shall further neglect 
mode coupling through parametric expiation in the an- 
alytic analysis, and explore the effect of mode coupling 
through numerical simulation. 

We therefore proceed by applying ((23)) with a single 
mode for some n into (|19p and then use Xn as a test 
function to extract the equation for ipn, i-e.: 

<i'rr'^Xn[{do(pn)Xn - Vn^Xn + VnfJ.'^Xn 

-eUoXniir) cos{Q.mt)ipnXn] = 0. (34) 
Using this yields 

d'^ipn 



dt^ 



+ UJn'fn - elloC Xnm COs{flmt)'^n = (35) 



where 



Xn 



= ^ I r XnXm 



Si(TO7r) ~ -Si(2n7r -I- mn) + -Si(2n7r — rmr). (36) 
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Modal equation including damping and parametric ex- 
citation — Incorporating both damping and parametric 
excitation, we arrive at the following equation 



dV„ 
dt2 



' dt 



+ ulifin - f-C^UoXnm COs(r2,„t)(/3„ 





(37) 



for each quasi-normal mode. To see the stability of these 
modes we cast ([37]) into the following canonical form for 
damped Mathieu equation: 



d^(Pn 
dT2 



2(^'^Vn_ ^ ^^^^ _ 2gcos(2T)(p„ = 
dr 



by using ([^0]) , ([^5]) , ([5^ , in terms of the following dimen- 
sionless quantities: 
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The stability domain near the principal parametric ex- 
citation frequency with a « 1 and q w have been ob- 
tained in [l7|. For a — 1, i.e. 0,„ = 2ti;„, (pn becomes 
unstable if the condition 



q 
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4 kI 



> 1 



(38) 



is satisfied. From we see that the above condition 
can be satisfied for sufficiently large (3. 

As an example, consider a model neutron star with 
radius R = 15km, equilibrium density po = 10^^ g/cm"^, 
which is radially pulsating with density wave speed v = 
0.75 c, amplitude parameter e = 1/3, mode index to = 12 
and frequency 0„/27r = 90 kHz. (Cf. (US]), ^ and 
([17]).) Using ([25)1 . this frequency is twice the frequency 
of the lowest quasi-normal mode of a massless scalar field 
with iiQ = 0, n = 1 and tOnjl-K = 45 kHz. It then follows 
from (|36p that Xnm = 7.5 x 10^*^, and so the estimated 
unstable (3 values using ([55]) yields (3 > 1400. 

Recall that the above estimated unstable (3 based on 
(155)1 does not take into account mode coupling, which 
could shift the boundary of instability moderately. To 
access this, we have performed numerical simulations via 
direct finite difference discretization of (fT9)l that allows 
all modes of (j) to be excited by the radial pulsation of 
the model neutron star. It is found that multiple mode 
excitations tend to lower the boundary of instability for f3 
by about one order of magnitude. Using the same data as 
the preceding example, unstable excitations of the scalar 
field are observed for /3 > 100 as show in Fig. [T] with 
13 = 100. 
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FIG. 1: Numerical simulation of the parametric excitation 
of scalar waves for f3 = 100 from inside a model pulsating 
neutron star where r < 1 and the resulting travelling scalar 
waves for r > 1. The vertical axis denotes r(j>{r,t) and the 
two horizontal axes denote radial distance r in unit of the star 
radius R and time t in unit of R/c. 
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